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A METHOD AND APPARATUS FOR DETECTING AND LOCATING NOISE 
SOURCES WHETHER CORRELATED OR NOT 

FIELD OF THE INVENTION 
The present invention relates to detecting and 
5 locating sources of noise in the general sense, using 

sensors that are appropriate for the nature of the noise 
source . 

The invention relates to a method of detecting and 
locating noise sources disposed in a space of one, two, 
10 or three dimensions and optionally correlated with one 
another, and presenting emission spectra of narrow or 
broad band. 

The invention finds particularly advantageous 
applications in the field of locating sources of noise 
15 optionally accompanied by echo and coming, for example, 
from vehicles, ships, aircraft, or firearms. 

BACKGROUND OF THE INVENTION 
In numerous applications, a need arises to be able 
to locate in relatively accurate manner a source of noise 
20 in order to take measures to neutralize it. Numerous 
solutions are known in the prior art for acoustically 
locating noise sources. The main known solutions make 
use of techniques for correlating signals delivered by 
detection sensors . 

2 5 Those techniques present the drawback of being 

particularly sensitive to interfering noise occurring in 
the environment of the measurement sensors. Furthermore, 
it must be considered that those techniques constitute 
specific methods that are adapted to each application 

3 0 under consideration. 

The technique in most widespread use involves 
antennas having a large number of sensors (several 
hundred) and a large computer system implementing beam 
forming so as to aim in a given direction in order to 
35 increase the signal-to-noise ratio. That method does not 
make any a priori assumption concerning the number of 
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sources and any possible correlation between them, which 
leads to a loss of resolution. 

OBJECTS AND SUMMARY OF THE INVENTION 
There therefore exists a need to have a general 
5 method of detecting and locating noise sources in space, 
when the number of noise sources is small and is known or 
overestimated. 

The invention seeks to satisfy this need by 
proposing a method of detecting and locating noise 
10 sources by means of sensors adapted to the nature of the 
noise source, the method presenting low implementation 
costs . 

To achieve this object, the method of the invention 
consists : 

15 - in taking the time-varying electrical signals 

delivered by the sensors, each signal s ± (t) delivered by a 
sensor being the sum of the signals Sj emitted by the 
noise sources; 

• in amplifying and filtering the time-varying 
20 electrical signals as taken ; 

• in digitizing the electrical signals; 

• in calculating the functional f, such that: 

det (<r k (co),r;(Q))> k, i = o to m) 

r (ni , • • • / z*j i • - - / n M) - — - — i 

det (< ^((oXTj (co)> k, 1 = 1 to M) 



with 



25 (T k (Q))) i =e 

<.,.> being the scalar product; 

. . c i being the vector constructed between the 
center of gravity of the sensors and the position of 
sensor i ; 

30 . . Tij being the unit vector in the direction 

defined by the center of gravity of the senors and source 

j: 

. . with T 0 = s; and 

. . with c = the speed of sound; and 
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• in minimizing the functional f relative to the 
vectors for j = 1 to M in such a manner as to determine 
the directions of the noise sources. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 Various other characteristics appear from the 

description given below with reference to the 
accompanying drawing which shows embodiments and 
implementations of the invention as non-limiting 
examples . 

10 Figure 1 is a diagram showing the principle of the 

detection method of the invention. 

Figure 2 is a diagram showing a detail 
characteristic to the method of the invention. 

MORE DETAILED DESCRIPTION 

15 As can be seen in Figure 1, the method of the 

invention consists in locating noise sources X lt X 2 , . . . , 
X jr . . . , X M where j varies over the range 1 to M, the 
sources being distributed in space and each emitting a 
respective signal Sj with j varying in the range 1 to M. 

2 0 The method of the invention consists in locating the 

noise sources Xj using sound wave or vibration sensors Y lf 
Y 2 , . . . , Yi, . . . , Y N where i varies over the range 1 to N, 
each delivering a respective time-varying electrical 
signal s x , s 2 , . . . , s ir . . . , s N . 

2 5 The method consists in taking the time -varying 

electrical signals s ± (t) delivered by each of the sensors 
and representative of the sums of the signals Sj emitted 
by the noise sources Xj . The signals s ± (t) received on 
the N senors on the basis of the sum of the contributions 

3 0 of the various sources is written as follows: 

s l(t) -. JU A ( t - at) 

where i = 1 to N, is the distance between the noise 

source X^ and the sensor Y if and c is the speed of sound 
in the ambient medium. 
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The term A ±j represents the attenuation due to 

propagation together with the sensitivity factor of the 

sensors and is expressed as follows: 

A ±j = B ± C(r ±i ) 

5 where i = 1 to N and j = 1 to M, where B i is the 

sensitivity coefficient of sensor Y i and where C(r^) is 
the attenuation coefficient due to propagation over a 
distance r ij7 - . 

The sensors Y ± are associated with respective 
10 electronic units (not shown) for amplifying and lowpass 
filtering the signals they pick up. The sensors are 
preferably matched in modulus and phase so that their 
sensitivities are identical- Thus, B i = G for i = 1' to N. 
Advantageously, in order to facilitate implementing 
15 the antenna of sensors as defined above, the sensors Y i 
are placed relatively close to one another. 
Consequently, for remote sources, the distance r ±j is of 
the order of the distance r^, i.e. the distance between 
the center of gravity of the sensors and the source . 
2 0 Thus, attenuation becomes a function of the distance 

only with CCr^) = C(r^) , with i = 1 to N and j = 1 to M. 
It can be deduced therefrom that : 



25 
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where i = 1 to N and j = 1 to M and: 



a ± (t) = ^(r^S^t - 



where i = 1 to N. 

Since the amplitudes of the sources X^ are unknown, 
the following equation can be written as follows, 
integrating the term a(r^) in S^i 




where i = 1 to N . 

Using Fourier transforms, the expression for the 
signals s ± (t) becomes: 
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(1) Si (CD) = S%.(co) .e c 

where i = 1 to N 

where s and S are the Fourier transforms of s and S 
respectively and where CO is angular frequency. 
5 This first equation (1) relates the received signals 

to the distance r ± ^ f i.e. to the positions of the sources 

As can be seen in Figure 2, other relationships can 
be expressed, associated with geometrical considerations 

10 enabling the distances to be related to the unit 

vector n jt which determines the direction defined by the 
center of gravity of the sensors and the source 
generating the signal S j . The position of the senors is 
defined by the vector C ± constructed from the positions of 

15 the sensors Y i and the position of their center of 

gravity. A development restricted to the first order of 
r ± j then provides : 

(2) r ±j ~ *j - <n j9 c ± > 

where i = 1 to N and j - 1 to M, and where <., . > is the 
2 0 scalar product. 

Thus, by replacing r ±j by the approximate expression 
given in (2) and integrating the phase term: 

-Jco— 

e 

which depends only on the source in the magnitude 
25 Sj(a)), equation (1) can be written: 

(3) s ± ((D) = X^(CO) .e 

where i = 1 to 

This relationship can also be expressed in matrix 
and vector form: 

30 (4) Si(CO) = £Sj(©) .Tj(©) 

j=i 



with, for ith coordinate of the vector T^z 
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-JO) ^ — — 



where i = 1 to W. 
Or indeed: 

(5) s(G>) = T.S(Q)) 
5 where T = matrix having the general term: 

-JO) — ^ 



T- ■ = e c 

In the presence of additive noise, equation (4) 
becomes : 

M 

( 6 ) ^ (CO) = X S j (CO) . Tj (CO) + B(CO) 

10 where B is the noise vector which depends on CO. 

The method of the invention consists in determining 
the directions of the sources defined by the vectors 
for j = 1 to M. 

When the sources are arbitrary, i.e. correlated or 
15 non-correlated, the probability of the presence of 
Gaussian noise B at the sensors Y ± is given by: 

where a and b depend on the variance. 

Thus, from equation (6) , the most probable position 
20 for the source is the position which minimizes the 
following magnitude: 

II 2 

s(to) = £.§_,■ (CO) .2V(C0) 

J=l II 

In other words, the projection of s onto the 
direction orthogonal to the hyperplane generated by the 
2 5 vectors Tj for j = 1 to M must be of minimum norm. 

That constitutes the square of the height of the 
parallelepiped constructed on the vectors s and iv, said 
height h being calculated as the ratio of the volume V to 
the base area S, i.e.: 

V 

30 h = — 

S 



The magnitudes V and S are expressed as a function 
of the determinants of the Gramm matrices in which the 
element (k, 1) is constituted by the scalar product: 

<T kf T* > 

with T 2 * being the vector that is the conjugate of T 2 . 
Thus , 

S 2 = det(<T JC# r 1 * > k,l = 1 to M) 

V 2 = det(<T^ / T I *> k, 1 = 0 to M) 
in which it is assumed T 0 = s. 
Thus : 




or indeed: 

f(n lf . . . , n jf . . . , n M ) = ||b|| 2 

This is a function of the vectors , which in three 
dimensions depends only on two angles 6^ and (p^ , i.e. 
elevation and azimuth, and on angular frequency CO. Any a 
priori knowledge about the spectra of the sources S j can 
also be used. For example: 

• for narrow band noise sources, provision is made 
to minimize the following functional f ± z 

(?) f x = X||Bto k )|| 2 

k 

co k being the angular frequencies of interest; and 

• for broad band noise sources, provision is made to 
minimize the following functional f 2 relating to the 
segments of interest : 

(8) f 2 = J||B(C0)fdC0 

In practice, instead of calculating 

Hf 

it is possible to use a sub-optimal method which consists 
in replacing the denominator S 2 by 1. It can then be 
shown, providing a search is being made for solution 
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directions that are sufficiently distinct, that 
solutions can be obtained that are close to those given 
by the exact method. The denominator, which cancels when 
at least two sources coincide, serves to eliminate 
5 interfering solutions where a plurality of directions are 
identical . 

When the sub-optimal method comprises broad band 
processing, it comprises minimizing the following 
functional f 3 : 

10 (9) f 3 = Jdet(<T^ / r*> k, 1 = 0 to M) dco 

This magnitude, which is expressed as a linear 
combination of the cross-correlation functions y i . of the 
signals s i and taken at points which are themselves a 
linear combination of the delays: 

15 

These cross-correlation functions are calculated 
only for delays having the same order of magnitude as the 
dimensions of the antenna divided by the speed of sound. 
Calculation can then advantageously be performed in the 

2 0 time domain as compared with calculation that is usually 

performed in the frequency domain on the basis of Fourier 
transforms . 

The above-described method thus consists: 

• starting from an antenna comprising a plurality of 
25 sensors (two to ten, and preferably two to five) in 

picking up acoustic or vibratory information; 

• in amplifying and in filtering the received 
signals so as to limit their spectrum and match the 
sensors in phase and in gain; 

3 0 • in digitizing the signals; and 

• in minimizing the projection of the vector s onto 
the direction orthogonal to the vectors Tj , using the 
following two calculation techniques a) and b) : 

a) + in obtaining the Fourier transforms of the 
35 signals s^- 
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+ in calculating one of the functionals f lf f 2 , 
using the above -de fined expressions for the Gramm matrix 
determinants; and 

+ minimizing one of the functionals f lt f 2 
5 depending on the number of sources to be located; 

b) + in using a simplified algorithm which 
consists in minimizing the functional f 3 ; 

+ in calculating the cross-correlation 
functions y±j * and 
10 + in minimizing a linear combination of the 

cross-correlation functions, depending on the number of 
sources . 

Once the minimization operation has been performed, 
the directions of the noise sources are determined. 
15 Advantageously, it is also possible to recover the 
characteristics of the noise sources . 

If N = M, i.e. if there are as many sensors as 
sources, then the system (5) can in general be inverted. 
If N > M, the problem can be reduced to a square 
20 system by premult iplying by: 

i.e. by the conjugate transposed matrix of T. System (5) 
then becomes : 

*T*.s(CD) = t T*.T.S(fi>) 

25 I.e. 

(10) S{C0) = ( t T*.T)- 1 . t T*.s(fl)) 

From equation (10), the signals Sj can be calculated 
so as to discover the characteristics of the sources . 
The description below gives an implementation for 
30 detecting one noise source (M - 1) using N sensors. 



This provides: 
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ll*MI 2 = 





N 


<n t Ci > 




+JC0 —±- 


||s(a))| 2 


I 


s ± * (co) e 




i=l 




N -J(o <D ' Ci> 






Xsi*(co)e 




N 


i=i 







N 



i.e. 



|B(CO)|| 2 = 



1_ 

N 



<n / c 1 -c k > 



(N - l ) .||s(o>)|| 2 - Xs^co) g^co) e DC °" c 

k,l 
k*l 



For broad band sources, it is thus a question of 
minimizing : 



j\\ B mf dco = 1 



(N-l)J ||s(co)|fdCD - XRe{J ^(CO) s^O)) e c " dco } 



3 co ± !? — 

c 



k*l 



I.e. writing the cross -correlations between the 
measured signals as y kl : 

(N - 1 )J ||s(C0)fdt - XYki ( c ) 

k,l 
k*l 

The problem thus reduces to maximizing: 

( <n. g, - ^> -j 

k*l 
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J||b(co)|| 2 dco = 



N 



by varying n which depends on only Np parameters, where 
Np is equal to the dimension of the space minus 1 unit. 

15 
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Without further elaboration, it is believed that one skilled in the art can, using the 
preceding description, utilize the present invention to its fullest extent. The preceding preferred 
specific embodiments are, therefore, to be construed as merely illustrative, and not limitative of 
the remainder of the disclosure in any way whatsoever. Also, any preceding examples can be 
repeated with similar success by substituting the generically or specifically described reactants 
and/or operating conditions of this invention for those used in such examples. 

Throughout the specification and claims, all temperatures are set forth uncorrected in 
degrees Celsius and, all parts and percentages are by weight, unless otherwise indicated. 

The entire disclosures of all applications, patents and publications, cited herein are 
incorporated by reference herein. 

From the foregoing description, one skilled in the art can easily ascertain the essential 
characteristics of this invention and, without departing from the spirit and scope thereof, can 
make various changes and modifications of the invention to adapt it to various usages and 
conditions. 



